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PROCEDE D'EXAMEN RADIOLOGIQUE D 1 UN OB JET 
DESCRIPTION 

Cette invention concerne un procede 
5 d'examen radiologique multi-energie d'un objet. 

Les procedes radiologiques consistent a 
faire traverser un objet a etudier par un rayonnement 
pour en deduire la repartition de differentes 
categories de materiaux, absorbant le rayonnement de 

10 fagons differentes, dans cet objet. Une application 
tres courante est 1 1 osteodensitometrie , ou on recherche 
la masse et la densite des tissus osseux dans un 
patient en distinguant ces tissus des tissus mous. 

II est usuel d'utiliser un spectre de 

15 rayonnement large et a le diviser en bandes mesurees 
separement par des canaux de mesure respectifs. Comme 
les coefficients d 1 absorption ou d 1 attenuation du 
rayonnement par l'une ou 1 1 autre des categories de 
tissu sont differents pour chacune des bandes, le 

20 probleme theorique se ramene a la solution d'un systeme 
d 1 equations dont le nombre est egal a celui des bandes 
de mesure et qui comprennent chacune deux inconnues 
(les epaisseurs ou les masses traversees de tissus mous 
et osseux) . La solution du probleme devient possible 

25 apres une calibration obtenue en faisant traverser par 
le rayonnement divers e talons dotes d ! epaisseurs 
connues de materiaux aux proprietes d' absorption 
analogues a celles des materiaux de mesure, et 
notamment le plexiglas et 1 1 hydroxyapatite pour simuler 

30 les tissus mous et les tissus osseux. On peut alors, 
calculer les parametres mathematiques d'un modele 
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reliant les mesures d 1 attenuation aux epaisseurs des 
materiaux. 

On notera qu'en realite le corps humain 
contient trois categories principales de tissus : les 
.5 tissus osseux, les tissus maigres et les tissus gras, 
mais qu'on n'en considere que deux en general a cause 
de difficultes a distinguer ces trois categories dans 
les mesures, de sorte qu'on confond volontairement les 
tissus maigres et gras. D'autres procedes sont ensuite 
10 appliques pour distinguer leurs proportions dans les 
tissus mous . 

La largeur du spectre permet de disposer 
d ! un nombre beaucoup plus grand de bandes de mesures 
qu'il ne serait necessaire pour obtenir une solution, 

15 et de les employer toutes pour obtenir des resultats 
plus precis en exploitant toute 1 1 information 
d f absorption obtenue . Dans 1 ! article "Measurement of 
bone mineral using a multiple-energy x-ray 
absorptiometry" , par J. Swanpalmer, R. Kullenberg, T. 

20 Hansson, Phys . Med . Biol , Vol 43, 1997 ;pp 379-387, ou on 
dispose de 23 bandes de mesures et ou on considere les 
trois categories de tissus, il est propose de combiner 
les groupes de mesures trois par trois de toutes les 
fagons possibles pour obtenir 1771 (23x22x21/2x3) 

25 systemes de trois equations a trois inconnues qui 
donnent autant de groupes de resultats. II convient 
alors d ! appliquer un critere de choix. Les auteurs 
conseillent de choisir comme resultat veridique celui 
qui est a la valeur moyenne, ou a la valeur mediane, 

30 pour le parametre le plus important, qui peut etre la 
masse osseuse traversee. 
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Contrairement a ce qu'on a conseille plus 
haut, les auteurs considerent les trois categories de 
tissus. Cela ne met pas en question la validite de leur 
methode, £ condition d'exposer le patient a une 
5 intensite d 1 irradiation beaucoup plus forte pour 
ramener les incertitudes des mesures a des proportions 
acceptables. II serait d'ailleurs possible de modifier 
simplement leur methode pour l ! appliquer a des mesures 
sur deux tissus, ou plus generalement deux categories 
10 de materiaux, seulement . Toutefois, il subsiste le 
defaut que cette methode ne tient pas compte du bruit 
sur les mesures et qu'elle conduit done a des resultats 
bruites . 

L 1 invention a pour objet un 

15 perfect ionnement de tels procedes par combinaison de 
resultats nombreux et comprend une amelioration du 
critere de choix des resultats. Plus precisement, elle 
concerne sous sa forme la plus generale un procede 
d'examen radiologique d'un objet oil sont considerees au 

2 0 moins deux categories de materiaux, comprenant ; 

l'emploi d'un rayonnement a spectre large ; des mesures 
du rayonnement par bandes du spectre ; des expressions 

(M) d'epaisseurs ou de masses des deux categories de 
materiaux traversees par le rayonnement, les 

25 expressions (M) etant des fonctions d f au moins deux 
des mesures (mesk) et de coefficient (A) ; et 
1 1 application d'un critere de choix parmi les 

A 

expressions ( M ) pour en deduire une expression 

(M finale) jugee veridique ; caracterise en ce que le 

3 0 critere de choix comprend un combinaison (f) des 

expressions avec des facteurs de ponderation (a) , et un 
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calcul de facteurs de ponderation de maniere que la 
combinaison ait un bruit minimal (variance minimale en 
langage mathematique) calcule d'apres les bruits sur 
les mesures (variances sur les mesures) . 
5 L 1 invention sera maintenant decrite en 

reference aux figures : 

- la figure 1 est une vue d ! un spectre, 

la figure 2 est une repartition de 

resultats , 

10 - et la figure 3 est un organigramme 

resumant 1 ' invention, auquel on pourra se reporter dans 
toute la description qui suit. 

L 1 attenuation du rayonnement peut etre 
exprimee par une fonction des epaisseurs traversees de 

15 chacun des materiaux d 1 indices x et y, ou de leurs 
masses M (densites par unite de surface) dans la 
direction du rayonnement. Le spectre de mesures de la 
figure 1 est divise en N bandes notees generalement par 
les indices i et j . Les attenuations varieront dans 

2 0 chacune des bandes en raison de coefficients 
d' absorption variables pour les deux materiaux. Si nous 
appelons mesi ou mesj les mesures pour une bande 
d'energie i ou j, les masses traversees M x et M y 
pourront chacune etre exprimees aussi, par exemple, par 

25 la formule generale 

M =Ai+A 2 .mesi+A 3 .meSj+A4 .mesf +A 5 .mesj +A 6 .mesi .mesj. . les 

mesures considerees dans cet exemple etant des meusures 
d 1 attenuation, on aura pour chaque canal de mesure i 
(correspondant a une bande du spectre) la relation mesi 

30 = In { no y .) ou noi estle nombre de phatons arrivant sur 
/ ni 

l'objet et ni celui des photons ayant traverse l ! objet. 
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Comme la non-linearite des fonctions M x et M y en 
fonction des mesures est faible en pratique, on peut 
s 1 accommoder de cette fonction polynomiale du deuxieme 
degre qui comprend six coefficients Ai a A € . Le degre du 
5 polynome peut etre ajustS en fonction du probleme ; Par 
exemple, pour l'analyse d'objets constitues de 
materiaux ayant des numeros atomiques plus eleves que 
celui des tissus biologiques, comme le controle non 
destructif pour 1 1 examen d'un objet metallique . les 
10 mesures considerees dans cet exemple etant des meusures 
d ! attenuation, on aura pour chaque canal de mesure i 
(correspondant a une bande du spectre) la relation mesi 



1' objet et ni celui des photons ayant traverse I 1 objet. 



une etape d'etalonnage a travers des etalons, appeles 
parfois fantomes ou cales dans l'art, et qui consistent 
en des pieces formees d'epaisseurs connues et 
differentes entre elles de materiaux simulant, par 

2 0 leurs proprietes d 1 attenuation, les materiaux de 
l'objet qui sera ef f ectivement ^ mesurer. Chacun des 
etalons est done soumis a une irradiation par le 
rayonnement pendant une longue duree qui permet de 
reduire 1 T influence du bruit sur les mesures. Le 

25 spectre mesure pour chacun des etalons donne encore N 
mesures resultant de la decomposition du spectre en 
autant de bandes. En combinant maintenant deux series 
de mesures i et j prises pour deux bandes et pour 
chacun des etalons, on recherche les coefficients A 

30 pour ajuster les fonctions M x et M y aux mesures. Dans le 
cas present, ou il y a six coefficients A pour chacune 




ou noi estle nombre de phatons arrivant sur 



15 



Ces coefficients peuvent etre trouves dans 
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des deux fonctions, et ou on considere les combinaisons 
de deux bandes de mesure, les mesures pour chacune des 
bandes devront, pour fournir une solution unique, 
porter sur six etalons . Des etalons plus nombreux 
5 pourront aussi etre utilises pour ameliorer la 
precision sur la solution. Une minimisation de fonction 
d'erreur sera alors appliquee. 

Cette determination des coefficients A est 
repetee pour diverses combinaisons de paires de 
10 mesures. II en etait encore ainsi dans 1 ! article 
anterieur mentionne plus haut ; cependant, on a 
constate qu'il etait inutile d'effectuer toutes les 

combinaisons, au nombre de pour exploiter 

completement les mesures et que (N-l) combinaison 
15 etaient en realite suffisantes pour recueillir toute 
1 1 information. 

Une fac^on preferee de proceder consiste a 
choisir au depart la bande de mesures qui est la moins 
bruitee (par exemple celle qui a le signal n - le 
20 nombre de photons a la reception- le plus important) et 
a l'associer successivement a chacune des autres bandes 
de mesures pour les combinaisons. On obtient finalement 
(N-l) estimations des deux fonctions M x et M y , qu'on 

note M Xi M 29 ...,MN -1 pour chacune de ces deux fonctions. 

2 5 A ce stade du procede, les fonctions M x et 

My, representatives des longueurs traversees de deux 
materiaux representatif s des tissus osseux et des 
tissus mous, pourront etre convert ies en fonction Mu, M v 
et M w representatives des longueurs equivalentes 

30 traversees des tissus osseux, des tissus maigres et des 
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tissus gras en combinant lineairement M x et M y de trois 
facpons differentes determinees par 1 1 experience . Ce 
procede de conversion est independant de l r invention, 
de j a connu dans 1 1 art • 
5 Comme il n'y a pas de raison de preferer 

I'une ou l 1 autre de ces estimations M, un critere de 
choix doit etre applique pour obtenir 1 1 estimation 

A 

M finale qui sera jugee veridique. Dans I 1 article 
anterieur, une des expressions obtenues etait 
10 directement select ionnee d ' apres un critere de 
classement (le resultat median) ou de moyenne des 
valeurs prises par les expressions pour un des 

resultats. Dans 1' invention, les expressions M seront 
combinees, par exemple lineairement, d ! apres la formule 

15 M finale = (al Ml) + (a2 M2)+...+ (a N-l M N -1) tout 
en minimisant le bruit ; les coefficients ai, etc. ont 
une somme egale a 1' unite (al + a2 + ....+ a N-l = 1) . 

Pour chaque canal de mesure, le bruit susr 
le nombre de photons suit une loi statistique de 

20 Poisson, dont le resultat est independant pour chacun 
des canaux. La matrice de covariance des N-l resultats 
peut etre exprimee d 1 apres la formule 

r = y dM, dM J 1 
iJ i~\dmes k dmes k N k 

La variance sur la combinaison lineaire 
2 5 donnant M finale s 1 exprime par la formule f = (ai , ... , a N - 
i) .r. fc (ai, a N -i) ; cette quantite f atteint une valeur 
optimale quand sa derivee selon toutes les variables 
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est nulle, c ' est-a-dire que 1 1 influence du bruit est 

r dflda x = 0 
df/da 2 =0 



minimisee, soit : « 



df/da^^O 



La resolution numerique de ce systeme 

5 fournit les coefficients al, a2 , etc et M finale, 
c 1 est-a-dire les masses traversees des deux categories 
de materiaux. 

Ce sont ces operations qui sont accomplies 
dans l f invention ; il est important de remarquer que 
10 les mesures (mes k ) employees dans la formule donnant T tj 

et ensuite sont les mesures effectuees a travers 
l ! objet a etudier lui-meme, mais pas celles qui ont ete 
faites a l'etalonnage pour determiner les coefficients 

A des fonctions M . 

15 L 1 invention permet d'echapper au dilemme de 

procedes anterieurs dans lesquels on renongait a une 
partie des energies de mesure ou au contraire on 
acceptait des mesures peu distinctes : elle utilise en 
effet tout le spectre, mais en le divisant en bandes 

2 0 assez nombreuses pour que la mesure de chacune puisse 
etre comparee de fagon utile a d'autres mesures, faites 
sur des bandes eloignees . Elle permet ainsi, entre 
autres avantages, d 1 etudier aussi bien des organismes 
maigres que gras . 

25 Enfin, l 1 invention peut etre generalisee a 

un nombre de materiaux plus grand que deux, ce qui peut 
etre interessant notamment pour les procedes d'imagerie 
a produit de contraste, ou trois variables doivent etre 
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considerees ou pour controle des bagages (recherche 
d'explosifs) / et elle peut etre appliquee en 
considerant des combinaisons autres que lineaires des 

fonctions M . 

5 Des application de 1 1 invention sont : 

. 1 1 osteodensitometrie 

- obtention de la densite de masse 

osseuse 

obtention de la composition 
10 corporelle (masse grasse, masse maigre) 

. le controle agro-alimentaire, par exemple 
la detection de fragments d'os dans de la viande ou la 
detection de morceaux de verres dans des plats 
15 cuisines, 

le controle des bagages : recherche 
d'explosifs, de produits illicites (armes, aliments, 
drogue...) 
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RE VEND I CAT I ONS 

1) Procede d'examen radiologique d'un objet 
ou sont considerees au moins deux categories de 

5 materiaux, comprenant : I'emploi d'un rayonnement a 
spectre large ; des mesures du rayonnement par bandes 

A 

du spectre ; des expressions ( M ) d'epaisseurs ou de 
masses des deux categories de materiaux traversees par 

le rayonnement, les expressions (M) etant des 
10 fonctions d ! au moins deux des mesures (mesk) et de 
coefficients (A) ; et 1 1 application d'un critere de 

A 

choix parmi les expressions (M) pour en deduire une 

A 

expression (Mfinale) jugee veridique ; caracterise en 
ce que le critere de choix comprend une combinaison (f) 
15 des expressions avec des facteurs de ponderation (a) , 
et un calcul des facteurs de ponderation de maniere que 
la combinaison ait un bruit minimal calcule d'apres un 
bruit sur les mesures. 

2) Procede selon la revendication 1, 
20 caracterise en ce que la combinaison (f) des 

expressions est lineaire. 

3) Procede selon la revendication 1, 
caracterise en ce que la variation de la combinaison 
est calculee avec une matrice de covariance ( IV ) des 

A 

25 expressions (M). 

4) Procede selon la revendication 1, 

caracterise en ce que les expressions (M) sont en 
nombre egal a celui des bandes, moins un, et etablies 
toujours avec une des bandes (io) et, respect ivement , 
3 0 chacune des autres bandes. 
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5) Procede selon la revendication 1, 
caracterise en ce que les coefficients (A) sont 
determines a une etape prealable de calibration. 

6) Procede selon la revendication 1, 
5 caracterise en ce qu ! il est applique a 

1 ' osteodensitometrie . 

7) Procede selon la revendication 1, 
caracterise en ce qu'il est applique a des controles 
agro-alimentaires . 

10 8) Procede selon la revendication 1, 

caracterise en ce qu ! il est applique au controle de 
bagages . 
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ABREGE DESCRIPTIF 



L 1 invention concerne un procede d'examen 
radiologique d'un objet ou sont considerees deux 
5 categories de materiaux, comprenant : 1 1 emploi d'un 
rayonnement a spectre large ; des mesures du 
rayonnement par bandes du spectre ; des expressions 

(M) d'epaisseurs ou de masses des deux categories de 
materiaux traversees par le rayonnement, les 

10 expressions (M) etant des fonctions d ! au moins deux 
des mesures (mes k ) et de coefficients (A) ; et 
1 1 application d'un critere de choix parmi les 

expressions (M) pour en deduire une expression 

(M finale) jugee veridique ; caracterise en ce que le 
15 critere de choix comprend une combinaison (f) des 
expressions avec des facteurs de ponderation (a) , et un 
calcul des facteurs de ponderation de maniere que la 
combinaison ait une variation minimale d'apres des 
variations des mesures. 

20 

Figure 3 . 
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